Transit times of radiation through photonic bandgap materials: Pulse reshaping of 
classical and quantum fields in the time domain 
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We study the propagation of model electromagnetic pulses through photonic bandgap materials 
and relate time-domain pulse reshaping to observable transit times. For layered dielectric mirrors 
we demonstrate how pulse reshaping of slowly varying classical fields results in transit-time delays 
that are equivalent to the group delay. The time-domain analysis emphasizes the causal origin of 
negative group delays. We also study an analogous fully-quantized model and show how the same 
time-domain analysis may be used to interpret observed delays of single-photon fields. 
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I. INTRODUCTION 



The propagation speed of particles through tunneling 
regions has been studied theoretically since the early days 
of quantum mechanics, and optical analogs of tunneling 
have been explored in recent years because of their ex- 
perimental accessibility. The physical meaning of group 
velocities greater than the vacuum speed of light c has 
been investigated in experimental work on the propaga- 
tion of photons through layered dielectric materials with 
photonic bandgaps at the frequency of the light [Q, |[ || . 
The measured delays of individual photons in these ex- 
periments were consistent with the group delay calcu- 
lated as the derivative of the phase of the transmission 
amplitude with respect to angular frequency. The group 
delay can be negative, seemingly consistent with propa- 
gation at a speed greater than c. As has been pointed out 
previously these anomalous delays are the result of pulse 
reshaping and they do not imply a violation of Einstein 
causality 0,|,|. 

In this paper we take advantage of the regularity of 
layered dielectric mirrors to explore explicit effects of 
pulse reshaping in the time domain on idealized model 
pulses, and we use this framework to interpret the origin 
of the observed pulse delays. We study reshaping of both 
classical field pulses and fully-quantized single-photon 
"pulses." We make a quantitative connection between 
the time-domain reshaping and the conventional group 
delay that yields insight into the origin of observed neg- 
ative delays. (Experimentally, the time-domain detec- 
tion of "super luminal" optical pulses is difficult, but such 
studies are possible in the terahertz M and microwave 0] 



regions of the electromagnetic spectrum, and analogous 
acoustic waves that travel faster than the speed of sound 
are experimentally accessible Q.) Our time-domain anal- 
ysis is a specific example of the "interference between 
. . . causally propagating consecutive components" that 
was proposed as universal mechanism for understanding 
delays in nondissipative media in Ref. , and we extend 
the classical ideas of Ref. Q into the domain of quantum 
mechanics and the transmission of single photons. Our 
analysis is complementary to analyses in the frequency 
domain like that of Ref. |^ . 

We consider propagation through mirrors comprised of 
layers of non-dispersive linear dielectric materials, each 
of which can be characterized by a single real index of re- 
fraction rii. This means that within a single material elec- 
tromagnetic waves propagate with equal phase and group 
velocities given by c/rii. The reduced speeds within such 
materials are partially responsible for observed pulse de- 
lays, but the reshaping due to multiple reflections also 
plays a significant role. Analogous anomalous sound 
speeds due to pulse reshaping in one-dimensional acous- 
tic bandgaps have also been investigated experimentally 

We study propagation of plane waves which are nor- 
mally incident on layered dielectrics; the thickness of a 
layer is given by di and the index by n^. To keep things 
simple we limit the analysis to materials in which the op- 
tical path length is the same for all layers, which means 
that the time it takes light to cross a layer is the same for 
all layers. For convenience we define a time Sa which is 
the round trip time within a layer, i.e., twice the transit 
time, 
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which is independent of the index i. We also define a 
time Sb which corresponds to the delay of a wavefront 
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propagating through the entire mirror: 



written 



Sb 



(2) 



In Sec.|l|we examine the propagation of classical fields 
and delineate the dependence of the observed delays on 
the two parameters Sa and Sb- In Sec. Ill we extend 
the analysis to fully-quantized single photon fields, and 
draw parallels between the propagation of classical and 
quantum fields. 



II. PROPAGATION OF CLASSICAL FIELDS 

Pulse reshaping in layered dielectrics is most dramatic 
for pulses with abrupt changes in amplitude, and we first 
consider the effects on square pulses. Although the re- 
shaping of such pulses leads to significant distortion, the 
simplicity of square pulses makes it easy to disentangle 
the effects leading to transit-time delays. After demon- 
strating the nature of the reshaping effects with square 
pulses we discuss slowly varying pulses, and derive quan- 
titative expressions for delays. 
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where in this expression (and throughout this article) A 
refers to the wavelength in vacuum. For a single dielectric 
layer the attenuation coefficients are 
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For more complex materials the coefficients Ti are con- 
structed from products of the appropriate single-interface 
reflection and transmission coefficients. 

For pulses with durations r that are large compared 
to the round-trip time Sa, many terms in Eq. (Q) have 
time to "turn on" before the pulse completely passes by 
the observation point. In the middle of such pulses the 
effective total transmission coefficient is 



N 



total 
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A. Reshaping of Square Pulses 

We consider a classical square pulse with complex mag- 
nitude £ and duration r, and assume that in the absence 
of any dielectric material this pulse arrives at the obser- 
vation point at t = 0, so that the observed complex field 
amplitude after propagating through a vacuum can be 
written 



for i < 
for < i < T 
for i > r 



(3) 



When the field passes through a layered dielectric slab 
the field arriving at the observation point is composed 
of multiple reflections. The leading edge of the trans- 
mitted field undergoes no reflections, and is delayed by 
the time Sb ; reflected fractions of the fleld will be addi- 
tionally delayed by multiples of Sa- We write the total 
field arriving at the observation point as the sum of terms 
grouped according the how much time the field has spent 
traversing the mirror. The first term corresponds to light 
that suffers no refiections; the second to light that spends 
an "extra" time Sa within the mirror; the third to light 
that spends an "extra" time 2Sa within the mirror, etc. 
(The terms are not grouped by the number of reflections. 
For example, some fields spending an "extra" time 25 
will have undergone two internal reflections, while some 
will have undergone four.) We separate the transmission 
coefficient for each term into a real factor describing the 
attenuation, and a complex factor giving the phase shift, 
so that the total field arriving at the observation point is 



where N is the number of terms that are "on" at the 
observation time. As N grows large the transmission ap- 
proaches that given by the standard steady-state trans- 
mission coefficient, which can be expressed as 
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Some examples of the effects of simple dielectric lay- 
ers on square pulses are illustrated in Fig. |l|. The top 
graph in the figure illustrates the time-dependence of the 
transmission through a single layer whose optical path 
length is (m-f 1/4)A, so that e*^'^'^''''/-^ = — 1, correspond- 
ing to a minimum in the steady-state transmission. The 
first arrival of the transmitted pulse is delayed by Sb, 
and this early arriving field is larger in magnitude than 
the rest of the pulse because at this time there are not 
yet any interference effects reducing the field. The lower 
graph displays transmission through a three-layer mir- 
ror. The additional layers result in an increased delay 
in the arrival of the leading edge and a reduction in the 
steady-state transmission coefficient. The additional lay- 
ers also result in an enhancement of the large magnitude 
of the leading-edge (when compared to the steady-state 
transmission). The relatively large early-time transmis- 
sion lasts for times on the order of Sa, the time between 
arrival of fields corresponding to terms in Eq. (jij). Note, 
though, that for multi-layer mirrors it takes more time 
for the effects of all of the multiple reflections to "turn 
on" at the observation point, and it takes longer for the 
transmitted pulse to settle down to the steady-state in- 
tensity. 
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FIG. 1: Transmission of square pulse through simple dielectric 
slabs. The top graph is for a single layer with n — 2.5 and 
the bottom is for three layers with m = 2.5, n2 — 1.25, 
and na — 2.5. In both graphs the optical path length of the 
dielectric layers is chosen to be (m + 1/4)A, where m is an 
integer, so that the steady-state transmission is a minimum. 
The time units correspond to the round-trip time through 
a single layer, i.e., Sa = 2nidi/c. The corresponding leading 
edge delays, Sb, are 0.3 and 0.7 in the upper and lower graphs 
respectively. 



The net effect of the layered dielectric is to create a 
transmitted field whose first arrival is delayed by 6b , but 
whose temporal intensity profile is reshaped so that there 
is relatively more intensity at early times compared to 
the intensity profile of the vacuum pulse. This shifting 
of intensity to relatively earlier times contributes to an 
effective advance of the pulse. If a single arrival time is 
to be assigned to a pulse (or a photon) , the relative role 
of these two effects must be accounted for when under- 
standing seemingly anomalous velocities. 

We note that for wavelengths corresponding to trans- 
mission maxima, the intensity profiles would differ from 
those illustrated in Fig. The time of first arrival would 
be unchanged, but the initial field would be relatively 
low, and would take time to "build up" to its steady-state 
value. The resulting pulse profile would show relatively 
more intensity at later times compared to the profile of 
the vacuum pulse, contributing to an effective delay. 



B. Transmission of slowly varying pulses 

For smoothly varying pulses any reshaping effects will 
be much smaller than those illustrated for square pulses 
in the previous section, but the transmitted pulse is built 
up from multiple reflections in much the same way. For 
pulses that vary slowly enough, the transmitted pulse 
that is constructed in this way will have the same shape 
as the incident pulse. We demonstrate that the transmit- 
ted pulse may be constructed such that it is delayed or 
advanced relative to the vacuum pulse. (No violation of 
causality is implied, and the intensity of the transmitted 
pulse is always lower than that of the vacuum pulse.) 

For simplicity we consider a portion of an incident 
pulse with linear amplitude modulation, and assume that 
the linear modulation has been in effect since a time to, 
so that the incident field can be written 



Ei{t)^£o[l + mit-to)]e- 



(8) 



In the following analysis we assume that the round trip 
time within the slab is much less than the time that the 
modulation has been in effect, i.e., 6a (t — to). This 
means that a very large number of the terms in a series 
like that of Eq. (|^) have "turned on," and the transmitted 
field is 
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where N is the number of terms that have "turned on." When N is large, the transmitted field is approximately 



l + m\ t — to — Sb— S/ 



(10) 



To first order in the small quantity 6a/ {t — to —5b) the transmitted field is 



t — to — Sb — SaRc 



(11) 



where a is a phase shift that will not be of further con- 
sequence in this analysis. 

Comparing the expression for the transmitted field 
given by Eq. ( |ll| ) to that of the incident field, Eq. (||), 
shows that the net result of the multiple reflections is an 
effective time delay of the linearly changing field given 
by 



■^effective 
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This effective delay can be positive, corresponding to 
a true delay, or it can be negative, corresponding to an 
advance. The exact value of the delay depends on the 
indices of refraction in the material comprising the mirror 
and the number of layers in the mirror. 

The first term in Eq. (|l2|), Sb, is simply the delay due 
to the reduced speed of wave-fronts within the dielectric 
materials. It is always a positive quantity, correspond- 
ing to an actual delay. The second term contains the 
more complicated effects of phased reflections, and may 
be positive or negative. If it is negative and greater in 
magnitude than Sb, the multiple reflection effects that 
led to the reshaping of square pulses dominate over the 
effect of reduced wave-front velocity, and the total de- 
lay is negative. It is important to note that the effective 
time delay does not arise from the simple shifting of the 
incident field at a given time to a new time. Rather, 
the effective delay is the result of the superposition of 
attenuated and phase-shifted fields from many previous 
times. 

Pulses will maintain their shapes and exhibit delays 
given by ^effective as long as the time scale characteriz- 
ing the modulation is long compared to Sa, the time be- 
tween the arrival of successive reflections. The peak of 
the transmitted pulse may arrive after the peak of the 
vacuum pulse would have arrived, or before. This is be- 
cause the transmitted pulse is not the result of simple 
attenuation of the incident peak, but rather it is con- 
structed from the superposition of many reflections, as 
in Eq. (^) . Slowly varying pulses are special in the sense 
that the newly constructed pulse has the same shape as 
the incident pulse. We emphasize that the effective delay 
does not apply to the arrival time of any feature associ- 
ated with an abrupt change in the field; the arrival of the 



leading edge of any disturbance associated with such an 
abrupt change will be delayed hy Sb- 

We conclude this section by demonstrating that the 
effective delay given by Eq. is identical to that pre- 
dicted by the conventional group delay, which is the 
derivative of the phase of the transmission amplitude 
with respect to angular frequency. Experimental mea- 
surements of delays have been consistent with the group 
delay, and our time-domain approach gives a physical 
picture of the origin of the observed delays. 

The steady-state transmission coefficient is given by 
Eq. d?]), and the phase of this transmission coefficient is 
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It is straightforward to show that 

del) 
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(13) 




Sb — SaRo 



which is equal to Jcffoctivc, the effective delay given by 
Eq. (|l^), which was derived above from time-domain con- 
siderations. 

The equivalence (for slowly varying pulses) of the 
group delay with the effective delay derived in the time- 
domain is a further demonstration that the group delay 
has a physical meaning, even in cases in which it results 
in a seemingly anomalous advance in the peak of the 
transmitted pulse relative to the peak of a vacuum pulse. 
For pulses which vary rapidly on the time-scale given by 
6a there will be significant distortion of the shape of the 
pulse that will depend on the details of the pulse shape 
and the characteristics of the dielectric mirror. Any dis- 
cussion of delays for such pulses must carefully account 
for such distortions in a way that is beyond the scope of 
the present analysis. 
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FIG. 2: Quantum model consisting of a pair of two-level 
atoms at fixed positions in a large multimode optical cav- 
ity. A dielectric mirror is centered in the cavity. We quantize 
the modes of the inhomogeneous cavity. The atom on the left 
begins in the excited state and is the source of a photon; the 
atom on the right serves as a detector. 



III. PROPAGATION OF SINGLE-PHOTON 
QUANTUM FIELDS 

In the preceding section we discussed the effects of lay- 
ered dielectric materials on the propagation of classical 
fields, but the experimental determinations of propaga- 
tion time cited above involved the detection of individual 
photons. In this section we demonstrate a way in which 
the time-domain picture presented for classical fields can 
be extended to single-photon quantum fields. We develop 
a model in which an excited atom spontaneously emits 
a quantized multimode photon, and we investigate the 
time-dependent probability for excitation of a detector 
atom located on the opposite side of a dielectric mir- 
ror from the emitting atom. The time dependent exci- 
tation probability displays interference effects that are 
exact analogs to those experienced by classical fields. 



A. Quantum Model 

We consider a large one-dimensional multimode opti- 
cal cavity of total length L which contains a symmetric 
dielectric mirror in the center of the cavity. (See Fig. |^.) 
The mirror is comprised of layers of homogeneous linear 
dielectric like that considered in Sec. |l[ We find the clas- 
sical standing-wave modes of the electromagnetic field in 
a cavity that includes the dielectric material, and quan- 
tize the modes of this inhomogeneous cavity. 

The cavity also contains a pair two-level atoms at fixed 
positions zi and Z2 on opposite sides of the dielectric re- 
gion. The atom at zi is initially in the excited state and 
spontaneously emits a photon into the quantized modes 
of the cavity; it serves as the source of the quantized field 
to be transmitted through the mirror. The atom at Z2 
begins in the ground state and serves as a detector of the 
transmitted radiation. (It is also possible to determine 
equivalent information about the transmitted radiation 
from quantities involving field operators such as the ex- 
pectation value of the intensity operator, (E~E~^) 



We use the excitation of a two-level atom because of the 
ease with which can compute this quantity to high pre- 
cision in our model.) 

The cavity is assumed to be large, in the sense that the 
length L is very much greater than the wavelength of the 
light emitted by the atoms, i.e., this is not a microcavity. 
The finite length L does not contribute to the physical 
phenomena under investigation; it simply provides a con- 
venient quantization volume for the field modes used in 
our calculations. 

We use a standard Hamiltonian of quantum optics to 
calculate the time evolution of the system, and pay par- 
ticular attention to the amplitude for the atom on the 
right side of the cavity to be found in the excited state. 
We note that the effects of the spontaneously emitted 
photon propagate causally in this model. The explicit 
form of the Hamiltonian we use is |ll|, |l^ 



interaction 



J2 hLo^^'^ai^^ + J2 ^^ma'^a,n 



+ XI ^ (5jmamcr+ ' + ffj>„aLc^- ' ) , (15) 



in which the atoms are labeled with the index j and the 
field modes with index to, and where 

t^(at) is 

the zero-field 

resonance frequency of both atoms, is the frequency 



of the 



,th 



field mode, 



and 



are the annihilation 



,tii 



mode, err', a 



U) Aj) 



and 



and creation operators for the to lii^j^K^, 

al are the pseudo-spin operators which act on atom j, 
and gjm gives the coupling of the j*"^ atom to the m}^ field 
mode. In this Hamiltonian we have made the standard 
electric-dipole and rotating-wave approximations. 

The mode frequencies LUm in the Hamiltonian are sim- 
ply those of the classical standing wave modes of the 
electromagnetic field. The spatial mode functions are 
normalized so that the energy per photon in the quan- 
tized modes is hujrn- The relative magnitudes of the cou- 
pling constants gjm reflect the spatial dependence of the 
classical mode functions, speciflcally the relative magni- 
tude of the mode functions at the positions of the two 
atoms. Calculation of the mode frequencies and spatial 
mode functions for the inhomogeneous cavity involves 
the solution of a classical boundary value problem. This 
is a straightforward process in principle, although the 
large number of boundaries in a multi-layer dielectric 
mirror leads to algebraic complexity. We used transfer 
matrix methods [ p3[ (adapted to electromagnetic stand- 
ing waves) and the "shooting method" 10 to determine 
numerically the mode frequencies. 

In the limit of a large cavity we may assume that the 
frequencies of all atomic transitions are very much greater 
than the fundamental frequency of the cavity. In this 
limit we can make the approximation that all modes that 
influence the dynamics of the system are near the atomic 
resonance, and the atom-field coupling constants can be 
factored into a product of a frequency-independent con- 



6 



stant and a space-dependent coupling factor. The cou- 
pling constants gj^ are given in terms of the electric 
dipole matrix element dj between the two levels of atom 
J, the effective volume of the cavity V, a mode-dependent 
normalization factor jV^i and the permittivity of free 
space eo, by 



gjm = ±djMm [ 2heoV J ^j )] 

= ±D,jJ\fm sin [krniL/2 - Zj)] , 



(16) 



where km is the wave- vector for mode rn, and in the last 
line we have defined the quantity 



(at) 



1/2 



2heoV 



(17) 



which is independent of the cavity mode-frequency. For 
symmetrically place atoms, modes with even spatial 
mode functions yield coupling constants with the same 
sign for each atom; odd mode functions give coupling 
constants of opposite signs. In performing our numerical 
calculations we use an equal number of modes above and 
below the atomic resonance frequency. 

The basis states for describing the system are 

• \e,g;0): left atom excited, right atom in ground 
state, no photon, 

• \g,e;0): right atom excited, left atom in ground 
state, no photon, 

• \g,g; Im)'- both atoms in ground state, one photon 
in m^^ cavity mode, 

and we write the state of the system as the linear com- 
bination 

mt)) = ci{t)\e,g;0} + C2{t)\g,e;0} 

+ Y,bm{t)\g,g;lm). (18) 

m 

In all the examples in this paper the system starts in the 
state 



IV-IO)) = |e,g;0). 



(19) 



and we pay particular attention to the complex ampli- 
tude C2{t) for the detector atom to be found in the ex- 
cited state. Although it is difficult physically to prepare 
a state which corresponds to our initial condition, this 
idealized state has the advantage that at t — Q all of the 
energy is localized at a single point (the position of the 
excited atom), making causal wavefronts evolving from 
this state particularly easy to identify. (A visualization of 
the propagation of the wavefronts of the intensity of the 
quantum field in similar models is presented in |l5l uM.) 



B. Method of Solution 

It is possible to find analytical solutions for the time 
evolution of atom-cavity systems with a single photon in 
simple inhomogeneous cavities ||l^, but the complexity 
of the mode structure for a cavity with a many-layered 
mirror makes this approach intractable. Therefore we 
construct numerical solutions for the coefficients ci(t), 
C2{t) and bk{t) of Eq. (|l|). 

We use the time-independent Schrodinger equation to 
determine the energies Eq and eigenstates \Eq) of the 
total Hamiltonian. The time evolution of the system is 
then straightforward to calculate. If the system begins 
in state 

\m) = |e,5;0) 

= 5]|£;,)(ii;,|e,g;0), (20) 

then the state of the system at a later time t is given by 

m)) ^Y.^-'''^"'''\E,){E,\e,g;^)- (21) 
q 

Projecting Eq. ( pi] ) onto the basis states gives the the 
time-dependent coefficients of Eq. (|8|): 

ci{t) = (e,5;0|V'(0) 



'iEqt/n 



[e,g-Q\Eq) 



(22) 



and 



C2{t) = (5,e;0|V'(i)) 

^ Y^e-^^^'l\g,e-MEq){Eq\e,g-Q), (23) 



hu{t) = {g,g-lkm)) 

^ Y^e-^^^'l\9,g-ME,){Eq\e,g:Q). (24) 



We use standard numerical matrix diagonalization rou- 
tines to determine the eigenvalues and eigenvectors used 
in these equations. We consider systems with as many 
as 2,000 modes, which leads to large matrix representa- 
tions of the Hamiltonian, but the the matrix is sparse. 
Our approach is similar to that used previously in several 
studies 0,111, 111 . 



C. Reshaping of single photon "pulses" 

The graphs of this section demonstrate the analogy 
between the classical field and the quantum amplitude 
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FIG. 3: Detector atom excitation probability in an empty 
cavity. The source atom decay rate is 71 = 16 and the detector 
atom rate is 72 = 256. The rapid rise at f = 1/2 is due 
to the arrival of the field at atom 2. The rate of the rise is 
determined by the detector response, and the slow exponential 
decay reflects the exponential shape of the pulse passing the 
detector atom. 



to find the detector atom in the excited state. The to- 
tal quantum amphtude is the result of interference from 
multiple reflections in very much the same way as the 
classical field is the result of multiple reflections. In this 
section we present results for three cavities: an empty, 
or vacuum cavity, a cavity with a simple homogeneous 
dielectric region, and a cavity containing a dielectric mir- 
ror. The results for the simple cavity help elucidate the 
more complicated behavior seen in the propagation of 
through the dielectric mirror. 

For atoms in an empty cavity with no dielectric mirror 
it is possible to find an analytic solution for the dynam- 
ics The probability for the source atom to be in 
the excited state decays exponentially with decay con- 
stant 71 = |riipL/c until the time at which reflections 
first interrupt the decay; for atomic positions zi — ~.2bL 
and Z2 = .25L this decay proceeds until reflections return 
to the atom time t — O.bL/c. For atoms at these posi- 
tions the amplitude to find the detector atom in the ex- 
cited state remains identically zero until time t = 0.5L/c, 
the time at which radiation can first reach the detector. 
Choosing units such that L/c = 1 (which we will use for 
the remainder of the article) the amplitude to find detec- 
tor atom in the excited state in an empty cavity is |T9] 



1\ 



71 - 72 



x(e-^(*-i)-e-*(*-4)), (25) 

in which the decay constant for the detector atom is 
72 = \ fl2\'^L/c, and the step- function expresses the turn- 



on of the excitation at t = 0.5 and the causal dynamics 
inherent in our model. 

In the limit of very fast detector atom response, i.e., 
72 ^ 00, the amplitude C2(t) will instantaneously re- 
fiect the strength of the field incident on the detector 
atom. Fig. ^ illustrates the detector atom response for 
the case 71 = 16 and 72 = 256. The rapid rise refiects 
the response of the detector atom to the sudden arrival at 
t = 0.5 of radiation from the source atom, and occurs on 
the time scale given by I/72. The slower decline reflects 
the exponential shape of the radiation pulse emitted by 
the source as it passes by the detector atom p^ . 
While in a given run of an actual experiment the excita- 
tion of the detector atom will be observed at a particular 
instant of time, the evolution of the excitation probabil- 
ity is continuous, and resembles the excitation expected 
for a classical oscillator driven by a classical field pulse. 

An illustration of the effect of a single dielectric slab 
on a pulse is given in Fig. ^. This figure gives the time- 
dependence of the detector atom excitation probability 
after transmission of the radiation through a single di- 
electric slab whose width and index are chosen so that 
the classical delays are Sa = 0.05 and Sb — 0.015 in the 
units of the figure. Because the decay rate of the source 
atom is relatively slow, this figure is analogous to the top 
graph in Fig. |l], which gives the intensity of a classical 
square pulse after passing through a single slab of dielec- 
tric material. The excitation of the detector atom "turns 
on" at the expected classical time, 0.5 -I- Sb, and is inter- 
rupted at multiples of Sa, the classical round-trip time 
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FIG. 4: Detector atom excitation probability in a cavity with 
a single homogeneous dielectric slab. The source atom decay 
rate is 71 = 1 and the detector atom rate is 72 = 1024. The 
width and index of the slab are chosen so that the classical 
delays are Sa = 0.05 and Sb ~ 0.015. The data points are 
the result of fully quantum-mechanical calculations, while the 
solid line is built from attenuated, delayed, and phase shifted 
versions of the vacuum response using classical delays, phase- 
shifts, and attenuation factors. 
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within the material. The excitation probabihty "settles 
down" to an attenuated and delayed version of the vac- 
uum "pulse" of Eq. (psj), where the attenuation factor is 
given by the classical iTg.s. 

The data points in the figure are the result of th e fully 
quantum mechanical calculations described in Sec. IIIB , 
while the solid line combines the quantum results for exci- 
tation in an empty cavity with the classical techniques for 
multiple reflections. In Sec. || the classical field is built 
up from a sum of appropriately attenuated, delayed, and 
phase shifted fields in Eq. (^); the solid line in Fig. ^ is 
built up from attenuated, delayed, and phase shifted ver- 
sions of the quantum excitation amplitude Eq. ( p5| ) using 
the same classical attenuation, delay, and phase shift pa- 
rameters that were used to produce the upper graph in 
Fig. I i.e., 

solid line < — > [c°(t- ^s)Toe*^'''=''«/^ 

-|-C^(i -Sb- 2(5A)r2e*27rc(5^+25B)/A 

+ •••]'. (26) 

The same principles apply to propagation through 
more complicated dielectric mirror structures. The ex- 
citation amplitude after transmission through an 11 el- 
ement mirror with alternating high and low index re- 
gions is illustrated in Fig. |^. This mirror is similar to 
the (HL)^H mirror used in the experiments of references 
[0, |[. The "high" index of refraction is nn — 3.0 and 
the "low" index is ul = 1.5, giving an intensity trans- 
mission coefficient at the minimum in the transmission 
of |Ts.s. P = 4.34 X 10"''. The upper graph in this figure 
also shows the empty-cavity excitation probability given 
by the square of Eq. ( p5[ ) scaled by the classical steady- 
state transmission factor, |Ts.s. p. 

The effects of the abrupt turn-on of the excitation are 
evident in the large variations in the early-time detection 
probability before the effects of multiple reflection have 
taken full effect. These variations are analogous to those 
in Figs. and ^ but they last longer in this case be- 
cause of the relatively large number of layers comprising 
the mirror. The lower graph in Fig. |^ displays the initial 
turn-on of the excitation of the detector atom. (Note that 
the vertical scale on this graph corresponds to amplitude 
rather than probability; the amplitude can assume nega- 
tive values just as the classical field E2 can.) As in Fig. ^, 
the data points are the result of the fully quantum me- 
chanical calculations, while the solid line combines the 
quantum results for excitation in an empty cavity from 
Eq. (^^ with the classical techniques for multiple reflec- 
tions as in Eq. (p6|). There is no excitation before the 
time t — 0.5 + Sb, and at later times the excitation is 
interrupted after successive multiples of the single-layer 
round trip time 6a- We call attention to the vastly differ- 
ent excitation probabilities at early and late times. The 
initial large peaks occur before the effects of multiple re- 
flection have reduced the transmission. 



The finite rise time evident in the lower graph in Fig. |^ 
is determined by the response of the detector atom char- 
acterized by 72. A more rapid detector response would 
show a more abrupt rise which more closely follows the 
step-function turn-on of the field at the detector atom. 
The modeling of more rapid changes would require the 
inclusion of more modes in our numerical analysis. 

The transmission of a photon from the source through 
a dielectric mirror and to the detector can happen via 
many indistinguishable pathways: it can travel directly 
without undergoing reflection; it may undergo a single 
reflection in one of many ways, or it may undergo mul- 
tiple reflections. The numerical results of this section 



3e-06 



^ 2e-06 
I 

X) 

p 

.2 le-06 



o 
X 




0.5 0.6 0.7 0. 

time (Units: L / c) 



0) 

3 



0.02 



0.015 



0.01 



0.005 



-0.005 



-0.01 




0.502 0.504 0.506 0.508 

time (Units: L/c) 



FIG. 5: Detector atom excitation in a cavity with an (HL)^H 
structure dielectric mirror. The source atom decay rate is 
ji — 4 and the detector atom rate is 72 = 1024. The width 
of the mirror elements is such that the classical delays are 
5a — 0.002 and Sb — 0.0057. The data points are the result 
of fully quantum- mechanical calculations, while the solid line 
is built from attenuated, delayed, and phase shifted versions 
of the vacuum response using classical delays, phase-shifts, 
and attenuation factors. See Eq. (0). 
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demonstrate that the excitation amphtude for the de- 
tector atom is built up from interfering ampHtudes for 
all of these processes in exactly the same way that the 
the classical transmitted field is built up from multiple 
reflections. 



IV. CONCLUSION 

We have examined the transmission of radiation 
through layered dielectric mirrors in the time domain. 
For slowly varying classical pulses we have derived a for- 
mula for pulse delays that takes into account the time- 
domain buildup of the steady-state transmission. The fi- 
nite time that it takes to build up this steady-state trans- 
mission results in pulse reshaping, and our formula helps 
delineate the competing effects of pulse reshaping and 
the reduced front velocities in dielectric materials. The 
delays calculated in our model are equivalent to those 
calculated from the conventional group delay, but our 
model provides an aid to understanding and interpreting 
the origin of the anomalous delay times that have been 
observed in experiments. 

We have also demonstrated that our interpretation can 



be extended to include quantized fields. We have studied 
numerically a model in which a spontaneously emitted 
photon propagates through a layered dielectric mirror, 
and excites a detector atom. The quantum amplitude to 
find the atom in the excited state is the result of inter- 
fering terms due to the possibility of multiple reflections 
in much the same way that the classical field is the re- 
sult of the interference of multiply reflected fields. The 
terms in the quantum amplitude "turn-on" at exactly the 
same time as the classical field terms, and with exactly 
the same relative amplitudes and phases. Individual pho- 
tons will be detected at a range of times described by this 
amplitude, and it is the distribution of arrival times of 
single photons that is shifted in time in exactly the same 
way as the classical pulse is shifted. 
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